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Linear maps of matrices describing evolution of density matrices for a quantum 
stem initially entangled with another are identified and found to be not always 
completely positive. They can even map a positive matrix to a matrix that is not 
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T^lj- , positive, unless we restrict the domain on which the map acts. Nevertheless, their 



form is similar to that of completely positive maps. Only some minus signs are 
inserted in the operator-sum representation. Each map is the difference of two 
completely positive maps. The maps are first obtained as maps of mean values and 



qh, then as maps of basis matrices. These forms also prove to be useful. An example for 



two entangled qubits is worked out in detail. Relation to earlier work is discussed. 



I. INTRODUCTION 

Linear maps of matrices can describe evolution of density matrices for a quantum system 
that interacts and entangles with another system 

HQS 

. The simplest case is when the 
density matrix for the initial state of the combined system is a product of density matrices 
for the individual systems. Then the evolution of the single system can be described by a 
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completely positive map. These maps have been extensively studied and used 
Here we consider the general case where the two systems may be entangled in the initial 
state. We ask what kind of map, if any, can describe the physics then. Completely positive 
maps can be used in quantum information processing because, with ability to decohere a 
system from its surroundings and initialize particular states, the two systems can be made 
separate so they have not been interacting and are not entangled when they are brought 
together in the initial state. What happens, though, when they are already entangled at the 
start? 

We find that evolution can generally be described by linear maps of matrices. They are 
not completely positive maps. They can even map a positive matrix to a matrix that is not 
positive. Nevertheless, basic forms of the maps are similar to those of completely positive 
maps. Only some minus signs are inserted in the operator-sum representation. Each map 
is the difference of two completely positive maps. These familiar forms follow simply from 
the fact that the map takes every Hermitian matrix to a Hermitian matrix. The maps are 
first obtained as maps of mean values and then as maps of basis matrices. These forms also 
prove to be useful. 

A new feature is that each map is made to be used for a particular set of states, to act in 
a particular domain. This is the set of states of the single system described by varying mean 
values of quantities for that system that are compatible with fixed mean values of other 
quantities for the combined system in describing an initial state of the combined system. 
We call that the compatibility domain. The map is defined for all matrices for the single 
system. In a domain that is larger than the compatibility domain, but still limited, every 
positive matrix is mapped to a positive matrix. We call that the positivity domain. We 
describe both domains for our example. 

To extract the map that describes the evolution of one system from the dynamics of 
the two combined systems, we calculate changes of mean values (expectation values) in the 
Heisenberg picture. This allows us to hold calculations to the minimum needed to find the 
changes in the quantities that describe the single system. To make clear what we are doing, 
we keep our focus on those quantities and keep them separate from the other quantities in 
the description of the combined system, which may be parameters in the map. 

There has been recognition of the limitations of completely positive maps in describing 
evolution of open quantum systems , but little effort has been made to use more general 
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maps there. Other considerations, including descriptions of entanglement and separability, 
have motivated substantial mathematical work on maps that are not completely positive 
but do take every positive matrix to a positive matrix ^,12, The maps we consider 
here do not need to have even that property. 

We begin with an example for two entangled qubits, which we work out in detail. Then 



we outline the extension to any system described by finite matrices. 
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n the concluding 



16| and point out 



section we discuss how what is done here relates to earlier work 
the errors in arguments that a map describing evolution of an open quantum system has to 
be completely positive. 

II. TWO-QUBIT EXAMPLES 

Consider two qubits described by two sets of Pauli matrices Si, S2, S3 and Hi, S2, H3. 
Let the Hamiltonian be 

H = ^ujE,E,. (2.1) 

The evolution of the S qubit is described by the mean values (Si), (S2) and (S3) at time 
zero changing to 

(e^-^*Sie-^^*) = (Si)coscjt- (S2Hi)sinu;t 
(e^-^*S2e~*^*) = (S2)coscjt+ (SiSi)sinu;t 

(e^^^Sse-^^*) = (S3) (2.2) 
at time t. These three mean values describe the state of the S qubit at time t. 

A. The basics for one time 

Look at this when utt is it/2. Then the mean values are changed to 

(Si)' = ai , (S2)' = a2 , (S3)' = (S3) (2.3) 

where 

ai = -(S2S1) , a2 = (SiSi). (2.4) 

We consider the oi, 02 to be parameters that describe the effect of the dynamics of the two 
qubits that drives the evolution of the S qubit, not quantities that are part of the description 
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of the initial state of the E qubit. What we do will apply to different initial states of the E 
qubit for the same fixed ai, 02- 

The change of mean values calculated in the Heisenberg picture determines the change 
of the density matrix in the Schrodinger picture. The density matrix 

p = i(l + (S)-S) (2.5) 

that describes the state of the E qubit at time zero is changed to the density matrix 

P = ^(1 + (E)' ■ E) = ^(1 + aiEi + asEs + (E3)E3) (2.6) 

that describes the state of the E qubit when ut is 7r/2. This is the same for all the different 
(E) that are compatible with the same fixed (E2H1) and (EiHi) in describing a possible 
initial state for the two qubits. We will refer to these as the compatible (E). 

To be meaningful, a map has to act on a substantial set of states. To insure that we 
have something substantial to consider here, we will assume that the set of compatible (E) 
is substantial. We will exclude those values of (E2H1) and (EiSi) that do not at least allow 
three-dimensional variation in the directions of the compatible (E). For example, we will 
not let (EiSi) be 1, because that would imply (E2) and (E3) are zero. The set of compatible 
(E) will be described more completely in Section FlI CI 

The change of density matrices can be extended to a linear map of all 2 x 2 matrices to 
2x2 matrices defined by 

1' = 1 + aiEi + a2E2 , E; = , E'2 = , E;^ = E3 . (2.7) 

This takes each density matrix p described by equation ()2.5|1 . for each compatible (E) in 
each different direction, to the density matrix 

p' = l(l' + (E)-S') (2.8) 

that is the same as that described by equation (|2.6p . This map takes every Hermitian matrix 
to a Hermitian matrix. It does not map every positive matrix to a positive matrix. 
The map takes 

^^=^(1 + ^3) (2.9) 

which is positive, to 

P' = ^(l + aiEi + a2E2 + E3) (2.10) 
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which is not positive. To see that P' is not positive, let 

ai = rcos9, a2=rsm9 (2-11) 

choose a vector ip such that 

{ip, lip) = 1 , {ip , T^i ip) = — r cos6'/vl + r^ 
{ip, = -rsinO/Vl + r^ , {ip , J:^ ^) = -1 / y/TT^ (2.12) 

and calculate 

iij, p» = i(l-v/rTH). (2.13) 

This is negative even when r is very small so that (S1H2) and (SiHi) are very small and 
there is room for a large set of compatible (S). 

Of course if p is a density matrix that gives a compatible mean value (S), the map takes 
p, described by equation ()2.5|1 . to the density matrix p' described by equation ()2.(i|l . which 
is positive. To see explicitly that p' is positive, consider that for any vector ip 

|(^, Si^)|2 + |(^, S2^)|2 + |(^, S3^)|2< |(^, (2.14) 

and if (S3) is compatible with (S1S2) and (EiSi) in describing a possible state for the two 
qubits, then 

{a,Y + (a2)' + (Ss)^ = (E2Si)^ + (SiSi)^ + {E^f < 1 (2.15) 

so that altogether 

\a,{^P, J:,^p) + a2{iJ, ^2^) + {^3)0^ , < , l^P). (2.16) 

The important difference between the density matrix p and the positive matrix P described 
by equation ()2.9|) is the factor (S3) multiplying S3 in the density matrix. If (S3) is changed 
to 1, the inequality ()2.15|1 can fail. The map can fail to take positive matrices to positive 
matrices when it extends beyond density matrices for compatible (S). 

The map can fail to be completely positive even within the limits of compatible (S) where 
it maps every positive matrix to a positive matrix. To see that, we extend the map to the 
two qubits by taking its product with the identity map of the matrices 1, Hi, S2, S3. We 
have used equations ()2.7jl to describe a map of 2 x 2 matrices. Now we use it to describe a 
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map of 4 X 4 matrices; each matrix in equations ()2.7|) is the product of the 2x2 matrix for 
the S qubit with the identity matrix for the H qubit. In addition we get 



(SiHfc)' = T.'{B!f. = 0, 
(E2Sfc)' = T^'^E'^ = 0, 

' [l-EkY =1'H; = (l + aiSi + a2S2)Sfc, (2.17) 



for = 1, 2, 3. This and the reinterpreted equations ()2.7|1 define a hnear map of 4x4 matrices 
to 4 X 4 matrices. If the map of 2 x 2 matrices defined by equations ()2.7|1 is completely 
positive, this map of 4 x 4 matrices should take every positive matrix to a positive matrix. 
We will see that it can fail to do that even when the 4x4 matrix being mapped is a density 
matrix for a possible initial state of the two qubits. 
If n is a density matrix for the two qubits then 



3 3 



n = ^ 1 1 + E(s.)s. + E(Sfc)2fc + E I (2.18) 

^ * j=l k=l j,k=l 



is mapped to 
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n' = M 1' + E(s.)s;- + E(s.)s'. + E (s,2.)(s,H,)' 1 . (2.19) 

V i=l k=l j,k=l J 

To test whether 11' is positive, let 

check that = to see that W is positive and is a density matrix, and calculate 

T..|n.., . i (: , ^ , 4 (i , <^ , , ,2,21) 

This holds if 11 is the density matrix for an initial state of the two qubits that gives the mean 
values — (S2H1) and (SiHi) used for ai and 02. We see that Tr[nW] can be negative. Both 
(SiHi) and (S3S3) are —1 for the state where the sum of the spins of the two qubits is zero. 
That state gives zero for (S), but nearby states will give an acceptable set of compatible 
(S) with Tr[nW] negative. 

The map is made to be used for the set of states, the set of density matrices, described by 
compatible (S) . We call that its compatibility domain. It includes all the initial states the E 
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qubit can have with the given (S2S1) and (SiSi). Outside the compatibihty domain, some 
density matrices are mapped to positive matrices, but others, including, for example, P from 
equation ()2.9j) . are not. Even inside its compatibility domain, the map is not completely 
positive. 

We can see that the compatibility domain is enough to give the linearity of the map 
physical meaning. Applied to density matrices, the linearity of the map says that if density 
matrices p and a are mapped to p' and a' then each density matrix 

T = qp+{l-q)(T (2.22) 

with < g < 1 is mapped to 

t' = qp' + (1 - q)a'. (2.23) 

Suppose p and a are density matrices for the S qubit that give mean values (S)p and (S)a-. 
If both (S)p and (S)o- are compatible with the same (S2S1) and (SiSi) in describing an 
initial state of the two qubits, then so is 

(S). = g(S)p + (l-g)(S).. (2.24) 

The compatibility domain is convex. Explicitly, if Hp and 11^ are density matrices for the 
two qubits written in the form of equation ()2.18|) with (S)p and (S)^- for (S) and the same 
(S2H1) and (SiHi), then 

= gnp + (i-g)n. (2.25) 

is a density matrix for the two qubits written in the same form with (S),- for (S) and the 
same (S2S1) and (SiHi). If p and a are in the compatibility domain, then so are all the r 
defined by equation ()2.22|1 . For these, the linearity described by equations ()2.22|1 and ()2.23|1 
has a meaningful physical interpretation. The compatibility domain will be described more 
completely in Section Hi CI 

A different map is an option if the initial state of the two qubits is a product state or if 
at least 

(S2H1) = (S2)(Hi) , (SiHi) = (2.26) 
Then the density matrix p' described by equation ()2.(ij) is 

P = \{l - (S2)(Si)Si + (Si)(Si)S2 + (S3)S3). (2.27) 
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This is obtained from equation ()2.8|) with the hnear map of 2 x 2 matrices defined either by 
equations ()2.7|) or by 

l' = l, S; = (Si)S2, S'2 = -(Si)Si, S'3 = S3. (2.28) 

With the latter, every positive matrix maps to a positive matrix. In fact the map is com- 
pletely positive. 

This completely positive map is defined by equations ()2.28j) for a given fixed value of (Hi). 
That puts no restrictions on (S), no limits on the initial state of the S qubit. Every (S) is 
compatible with any (Si) in describing an initial state of the two qubits for which equations 
fl2.26|) hold; every state of the S qubit can be combined with any state of the S qubit in a 
product state for the two qubits. However we will see that the (S) compatible with given 
nonzero (S22i) and (SiSi) in product states for the two qubits fill only a two-dimensional 
set embedded in the three-dimensional compatibility domain. 

The completely positive map defined by equations ()2.28|1 is an option only when equations 
fl2.26|) hold. Then both maps, from equations ()2.7|) and ()2.28|) . reproduce the evolution of 
the S qubit. There is a map defined by equations ()2.7p for almost every initial state of the 
two qubits, with (S2S1) and (SiSi) changing continuously from state to state. Switching to 
the completely positive map when equations ()2.2(i|l hold would be a discontinuous change. 

B. Time dependence 

From the mean values ()2.2|1 for any t, the same steps as before with equations ()2.5|1 . ()2.fj|l 
and (123) yield 

1' = 1 + (aiSi + 02^2) sin cut 
= J^icosut , S2 = S2 cos cut , S3 = S3. (2.29) 

This defines a linear map Q ^ of all 2 x 2 matrices to 2 x 2 matrices described by 

Q'rs = ^^rj;skQjk (2.30) 

jk 



with 



1 



a sin cut 








^a* sin out cos out ^ 








V 



cos out sin out 



■^a* sin out 

1 



where a = ai + ia2 and the rows and columns of B are in the order 11, 12, 21, 22. 
A vector 

^ A 



1 or 3 — 



^a* sin cut 
I a sin out 
X 



is an eigenvector of B with eigenvalue A if 



A H — lap sin^ out + X cosa;t = A^. 
4 



This yields two eigenvalues 



Ai = - + cos ujt + y (1 + cos ujty + {al"^ sin^ ut 
A3 = ^ + cosa;^ — + cosa;^)^ + |aPsin^a;t^ 

and eigenvectors 

"01 or 3 = i>i for A = Ai 
= ip3 for A = A3. 

Note that V'l and ■03 are orthogonal because 

A1A3 = —-\af sin^ out. 

The squares of the lengths of the eigenvectors are 

2{Xi + \\a\'^ sin'^ out 



for n = 1, 3. A vector 



= 2An(l + cosa;i) + |a|^ sin^a;f 

( X ^ 

— ha* sincL^t 

"02 or 4 = 



V 



I a sin ojt 
-A 
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is an eigenvector of B with eigenvalue A if 

1 



X + ^\a\'^sm^ujt- Xcosujt = X'^. (2.39) 



This yields two eigenvalues 



1 



A2 = - ( 1 — coscut + Y (1 — coscut)^ + |ap sin^ cut 
1 



and eigenvectors 



A4 = ^ ( 1 -coswt- -costut)2+ laPsin^cut) (2.40) 



■^2 or 4 = 'ip2 for A = A2 

= for A = A4. (2.41) 



Note that ip2 and V'4 are orthogonal because 

1 



A2A4 = -^lal^sin^cjt. (2.42) 

The squares of the lengths of the eigenvectors are 

ll^nlP = 2(^A2 + l|apsin2cut) 

= 2An(l - coscut) + lapsin^tjt (2.43) 

for n = 2, 4. 

We see that, in all but a few exceptional cases, B has two positive eigenvalues Ai and A2 
and two negative eigenvalues A3 and A4. That means the map is not completely positive; for 
a completely positive map, i3 is a positive matrix and its eigenvalues are all non-negative. 
A plot of the eigenvalues of i3 as a function of ut when | a p is 1/2 is shown in Figure ^ The 
two negative eigenvalues A3 and A4 go to zero when ut is nvr; the map is the identity map 
for even n and rotation by vr around the z axis for odd n. 

The spectral decomposition 

4 

B=J2Xr^\n){n\ (2.44) 



n=l 



with 



yields 



l'^) = ^iV'n) (2.45) 



Brj-sk = K{rj\n){sk\n)* = ^ sign(A„)C(n)rjC(?2)),^ (2.46) 



n=l n=l 
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FIG. 1: (color online) The eigenvalues of i3 as a function of ujt when |ap is ^. The dot-dash (red) 
line is Ai, the solid (green) line is A2, the dashed (blue) line is A3, and dotted (black) line is A4. 

with 



C{n)rj = ,/\K\{rj\n) = ttt-^^^-^'I^") 



so equation ()2.30|) is 



ll^nl 



Q'rs = sign(A„) ^C(n)^jgjfcC(n)l, 

n=l jk 



or 



Q' = ^sign(AJC(n)QC(n)t. 

n=l 

Since TrQ' = TrQ for all Q for our map, 

^ sign(A„)C(n)tC(n) = 1. 



(2.47) 



(2.48) 



(2.49) 



(2.50) 



n=l 



Except for the minus signs, these equations are the same as for completely positive maps. 
Explicitly we have 



C{n) 







2A„(1 + cos cut) + 


a 


^ sin^ ujt 



Xn + ^(ai^i + 02^2) sin cut 



(2.51) 



for n = 1, 3, and 



|A„ 




2A„(1 — costut) + 


a 


^ sin^ ut 



C{n) 



for n = 2,4. For small ut and nonzero \a\ 



A„S3 + -(02^1 - 01^2) sin cut 



(2.52) 



Ai = 2 - hutf + l\a\\u;tf 

/ o 
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and 



C(l) 
C{2) 
C{3) 
C{4) 



A2 

A3 
A4 



1, , 1/ n2 1 



1-^ + ^(01^1 + 02^2) 




2 a 



S3 + 



8|a| 



(a;i)2(ia2Si - iaiE2) 



|ar . .^2 



+ — (oiEi + 0222) 




S3 + 



8|a| 



(2.53) 



(a;t)2(ia2Si -^1^2). (2.54) 



C. Compatibility and positivity domains 



Now we describe the compatibility and positivity domains completely and precisely. To 
write equations for the compatibility domain, we make a convenient choice of components 

— * 

for (E). Suppose oi and 02 are given. Then (EiSi) and (S2S1) are fixed. Let 

(SiSi)Si + (S2Si)S2 



(SiSi)2+(S2Si)2 

(S2^i)Si — (Si^i)E2 



^(EiSi)2 + (S2Si)2 • 
Then E+ and E_ anticommute, their squares are both 1, and (E_Si) is zero, 



and 



(2.55) 



(2.56) 



(2.57) 



(Si'^i)Si'^i + (S2i^i)E2.::ii — (E_|_,=ii)E_|_,:ii. 

The compatibility domain is the set of (E), or (E+), (E_), (E3), that are compatible with 
the given (E+Si) and zero (E_Si) in describing a possible initial state for the two qubits. 

Basic outlines of the compatibihty domain are easy to see. When (E_)_) is zero, the 
compatibility domain includes the (E_), (E3) such that 



(E_)^ + (E3)^ + (E+Si)^<l 



(2.58) 
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because for these 

n = ^(1 + + (S3)S3 + (2.59) 

is a density matrix for the two qubits. Larger and (S3) are not included. If 

(x_)2 + (x3)^ + (S+Si)2 = l (2.60) 



and r > 1, then 



n = ^(l + rx_S_+rx3S3 + (S+Si)S+Si 

3 33 
+ Vj-i + ^3iS3Si + ^jk^j^k) (2.61) 

j=l j=l k=2 



is not a density matrix for any yj and zjk because 



W=^il- - X3S3 - (S+Hi)S+Hi) (2.62) 

is a density matrix and 

Tr[nW] = ^(1 - rix.f - r(x3)2 - (S+Si)^) < 0. (2.63) 

When is zero, the compatibihty domain is just the circular area described by (12.58^ : it 
cannot be extended in any direction described by any ratio of (S_) and (S3). This projection 
of the compatibility domain on the (S_), (S3) plane is shown in Figure 121-A for the case 
where (S+Si) is 

When (S3) is zero, the compatibility domain is the elliptical area of (S_), (S+) such that 

(S_)^ 



1 - (S+Si)2 

To see this, we find when the eigenvalues of 



+ (S+)' < 1. (2.64) 



n = ^(1 + (S_)S_ + (S3)S3 + (S+)S+ + (S+Si)S+Si + (Hi)Si + (S3Hi)S3Si) (2.65) 
are all nonnegative so that 11 is a density matrix for the two qubits. Let 

n= ^(l + (Si)Si + M). (2.66) 

Then 

M' = (S_)^ + (S3)^ + (S+)2 + (S+Si)2 + (S3Si)2 + 2(S3)(S3Hi)Hi + 2(S+)(S+Si)S,. (2.67) 
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The eigenvalues of M are the square roots of the eigenvalues of M^. When Si has eigenvalue 
+ 1, the eigenvalues of 11 are 

\{l + {^,) ± ^m^)) (2.68) 

where m^(+) is with Si replaced by its eigenvalue +1. When Si has eigenvalue —1, the 
eigenvalues of 11 are 

^(l-(2i)±^A^) (2.69) 

where m^(— ) is with Si replaced by its eigenvalue —1. The eigenvalues of 11 are all 
nonnegative if 

m'{+)<{l + {E,)r (2.70) 

m2(-)<(l-(Si))2 (2.71) 

and (Si)^ < 1. When (S3) is zero, the areas of allowed by the inequalities (j2.7(jp 

and ()2.7H) are largest when (S3S1) is zero. Then as (Si) varies from —1 to 1 the inequalities 
(I2.70|) and ()2.71|) describe the area of an ellipse with foci at ±(S_|_Si) on the (S_|_) axis; they 
say that the distance from a point with coordinates (S_), (S+) to the focus at — (S+Si) is 
bounded by 1 + (Si) and the distance to the focus at (S+Si) is bounded by 1 — (Si), so 
the sum of the distances is bounded by 2. That gives the elliptical described by fll^ . 
We conclude that it is the compatibility domain when (S3) is zero. This conclusion is not 
changed if 11 is given additional terms involving S2, S3, SjS2, SjS3. Each eigenvalue that 
we considered is a diagonal matrix element [ip, Hip) with ip and eigenvector of Si as well 
as an eigenvector of the 11 we considered, so (■?/', S2-?/'), (■?/', S3'?/^), {ip,T,jE2ip), (i^jT^jE^ip) are 
zero. Additional terms will change the eigenvalues and eigenvectors of 11 but will not change 
the diagonal matrix elements we considered. They have to be nonnegative if 11 is a density 
matrix. That is all we need to show that the inequality ()2.fj4j) describes the compatibility 
domain when (S3) is zero. The projection of the compatibility domain on the (S_|_), (S„) 
plane is shown in Figure 121-B for the case where (S+Si) is l/\/2. 

When ai and 02 are not both zero, all the product states for the two qubits that are 
compatible with the given (S+Si) and zero (S_Si) are for (S) in the projection of the 
compatibility domain in the (S3), (S+) plane. If 

(S_)(Si) = (S_Si) = 

(S+)(Si) = (S+Si)^0 (2.72) 
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then = and 

> (S+Hi)^ (2.73) 
There is a compatible product state for each such (S4.) and each (S3) such that 

(Sa)' < 1 - (S+)2, (2.74) 

with (S_) = 0. The (S) for compatible product states fill the two areas in the (S+), (S3) 
plane bounded by sections of the unit circle from ()2.74j) and straight lines from ()2.73|1 . These 
areas are shown in Figure 121-C for the case where (S+Hi) is l/\/2. 

Since (S) cannot be outside the unit circle for any state, these sections of the unit circle 
are on the boundary of the compatibility domain. We can conclude that the boundary of 
the projection of the compatibility domain in the (S_|_), (S3) plane is completed by straight 
lines with constant values of (S3) between the sections of the unit circle, because we proved 
the compatibility domain is convex and from (j2.58j) . ()2.74|1 and ()2.73j) we see that (S3)^ 
cannot be larger when (S4.) is zero then it is at the termini of the sections of the unit circle. 
The complete boundary in shown in Figure 121-C for the case where (S+Si) is l/-\/2- 




FIG. 2: Sections of the compatibility domain when (S+Hi) = The area enclosed by the thick 
solid line is the compatibilty domain. The dotted line shows the unit circle. The shaded area in 
(C) shows the (S) for product states compatible with the given (S+Hi) and zero (S_Hi). 

We will show that the compatibility domain is the set of (S) where 

^((S_)2 + (S+)2 + (S+Si)2)2 - 4(S+)2(S+Si)2 < 2 - 2(S3)^ - (S.)^ - (S_,)2 - (S+Si)l 

(2.75) 
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First let us see what this says. Squaring both sides of ()2.75|) gives 



{^-r + + (S3)^ + (S+SO^ - < 1- (2.76) 

When is zero, ()2.76p is the inequahty ()2.58|) that describes the circular projection of 
the compatibility domain in the (S3) plane. When (S3) is zero, ()2.76|) is the inequality 

()2.64|1 that describes the elliptical projection of the compatibility domain in the (S_), (S+) 
plane. If (S3)^ is between zero and 1 — (S+Si)^, then ()2.7(')j) is 

+ < 1- (2.77) 



1-(S+Si)2-(S3)2 1-(S3)2 

A contour of the compatibility domain at constant (S3) is an ellipse. As (S3)^ approaches 
1 — (S+Hi)^ the semi-minor axis shrinks to zero and the semi-major axis goes to (S+Si), 
so the ellipse reduces to a line from — (S+Si) to (S+Hi) along the (S+) axis. When (S_) is 
zero, ()2.75p is 

(S3P<l_<5±)l±<5±£l>!_K5±>!^(5±£in (2.78) 
which is ()2.74|1 when (S^)^ > (S+Si)^ and is 

(S3)' < 1 - (S+Si)2 (2.79) 

when (S_,_)^ < (S+Si)^. That describes the area bounded by sections of the unit circle and 
straight lines that is the projection of the compatibility domain in the (S.,.), (S3) plane. 

When (S+Si) is zero, ()2.75|1 just says that (S) is on or inside the unit sphere; then there 
is no restriction on (S) from compatibility. A three-dimensional view of the compatibility 
domain is shown in Figure Elfor the case where (S+Si) is 

The inequality ()2.75p puts a bound on (S3)^ for each (S_) and (S+). In particular, it 
says (S3)^ can never be larger than the values it has when (S_) is zero; the bound ()2.79|l 
holds for the entire compatibility domain. For (S3)^ within this bound, the left side of ()2.7(j|l 
is an increasing function of (S+)^. The inequality ()2.7fi|) puts a bound on (S_)^ for each 
(S+) and (S3) and a bound on (S+)^ for each (S_) and (S3). 

To show that the set of (S) described by (I2.75P is in the compatibility domain, we show 
that for each (S) that satisfies ()2.75|) there is a 11 described by ()2.65|) that is a density matrix 
for the two qubits. We let 
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FIG. 3: (color online) The compatibility domain generated using Mathematica for the case where 
(E2H1) and (SiHi) are both i. The dotted sphere is the unit sphere (the Bloch sphere) that 
represents all possible states of the qubit. 

and 

(S3Si) = (S3)(Hi). (2.81) 
Then the inequahties dTTHll and (ETTD are both (IT7T)|l . From (IT7T?|l . which (ETHll imphes, 

l(2i)|<J^^<l (2.82) 
(^+-1/ 

for < and from (071) . which holds for any (S), 

I(H,)I<^<1 (2.83) 

for > (S+Hi)^. This implies that the eigenvalues of 11 are all nonnegative, which 

means 11 is a density matrix for the two qubits. 

The inequality (j2.76|) by itself does not imply that (S) is in the compatibility domain. 
The equality limit of (j2.76|) is a quadratic equation for (Ss)^. The equality limit of (|2.75p 
is one solution. In the other solution, the sign of the square root in ()2.75p is changed. That 
changes the sign of the term with the absolute value in ()2.78|) , which extends the boundary 
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to include the entire area of the unit circle in the (S3) plane. The bounds ()2.79|) on 

(£3)^ and ()2.82j) on | (Hi) | do not hold for the other solution. They are not implied by ()2.75p . 

We have shown that the set of (S) described by the inequality ()2.75j) is in the compatibility 
domain. The compatibility domain is the same for all t. In a larger domain, which we call 
the positivity domain, every positive matrix is mapped to a positive matrix. The positivity 
domain depends on the time t. We will show that the set of (S) described by the inequality 
fl2.75|) is also the intersection of all the positivity domains for different t. That implies it is 
the compatibility domain; the compatibility domain cannot be larger, because it must be in 
every positivity domain for every t. 

The positivity domain for each t is easily found from the map of mean values 

(51) ' = (Tii) cos out + ai sin Lut 

(52) ' = (S2) cosa;t + 02 sina;t 

(53) ' = (S3). (2.84) 

Regardless of whether (S) is compatible, the density matrix for (S), described by equation 
(12. 5|) . is mapped to a positive matrix, which is the density matrix for (S)' described by the 
first half of equation ()2.6p . if 

((Si)')' + ((S2)')' + ((S3)0'<1 (2.85) 

— * 

which means (S)' is on or inside the unit sphere described by 

(Si)' = sin^ cos(/9 , (S2)' = sin6' sin , (S3)' = cos^ (2.86) 

with 6, ip varying over all directions. Then (S) is on or inside the surface described by 

,^ , sin 6* cos 09 

(51) = — aitancjtH 

cos cut 

sin ^ sin 

(52) = —a2ta,nujt-\ 

cos cut 

(53) = cose (2.87) 

which is obtained from the unit sphere by moving the center distances —ai tan tut and 
—02 tanut in the x and y directions and stretching the x and the y dimensions by a factor of 
1/ cos cut. The positivity domain is the intersection of this surface and its interior with the 
unit sphere and its interior, since (S) must also be on or inside the unit sphere. The posi- 
tivity domain for different values oi ut is shown in Figured When cut is 7r/2, the restriction 
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(A) (B) (C) (D) (E) 




FIG. 4: (color online) The positivity domains for (left to right) = %^ % and ^ when 
oi is — ^ and 02 is ^. The surface of the unit sphere is shown with dotted lines where it is not the 
surface of the positivity domain. When Lot is 0, the positivity domain is just the whole unit sphere. 

(ITHHll is just that 

(S3)2 < 1 - {a^f - {a^f. (2.88) 

Then the positivity domain is the part of the unit sphere where (Ss)^ is within this bound. 
If Oi and 02 are not both zero, and t is not zero, the positivity domain does not include the 
north pole point that corresponds to the matrix P of equation ()2.9p . 

If Oi and 02 are both zero, the positivity domain is the entire interior and surface of 
the unit sphere. Then the map takes every density matrix to a density matrix and every 
positive matrix to a positive matrix. In fact the map is completely positive for all t. The 
two eigenvalues of B that are generally negative, A3 and A4, are zero, so C(3) and C(4) are 
zero. That leaves two positive eigenvalues 

Ai = l + cosu;t, A2 = l — coscut (2.89) 

and just 

= ^l±Spl, m = /l^S3. (2^90) 

Consider three sets: the intersection of all the positivity domains for different t, the 
compatibihty domain, and the set of (S) described by the inequality ()2.75|1 . We know these 
sets are nested; the intersection of the positivity domains contains the compatibility domain 
because every positivity domain contains the compatibility domain, and we showed that the 
compatibility domain contains the set of (S) described by ()2.75|) . Now we will show that 
these three sets are the same; we will show that every point on the boundary of the set of 
(S) described by ()2.75|) is also on the boundary of a positivity domain for some t. 
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In terms of the used to describe the compatibihty domain, the equations 

()2.87|) for (S) on the boundary of the positivity domain for time t are 

/V \ ^^^^ • f ^ 
S+ = -sm{ip-a) 

cos ivt 

sin 6 

= (S+Hi)tana;t cosfoj — a) 

cosut 

(S3) = cos 9 (2.91) 



with 



If 



then 



ai = (S_|_Si) cosa , 02 = (S+Si) sin a. (2.92) 



sin ut = = {^^^i)cos{v-a) ^2.93) 



(S+) = -sin^sin/3 = - - {Es)^smp 

(S_) = -v/sin2e-(S+Si)2cos/3 = -^1 - (Sg)^ - (S+Hi)2 cos/5 (2.94) 



. sin((/?-a) 
smp = 



cos cut 



where 



/sin - (S+Si)2 cos(¥? - a) 

cos/3 = — 

sm u cos uot 

sin 6 

tan/7 = =tan(y9 — g). (2.95) 

V'sin2^-(S+Si)2 

You can check that the sum of the squares of the formulas for sin/5 and cos (3 is 1, so the 
designations sin/5 and cos/5 are allowed. Each (E) described by these equations is on the 
boundary of a positivity domain. Equations (j2.94j) also describe the ellipses of (j2.77p that 
are the contours of the boundary of the set of (S) described by the inequality ()2.75|) . From 
equations ()2.95p we see that all values of (3 from to 27r are included as ip — a varies from 
to 27r, so the whole of each ellipse is included. The bound ()2.79|) on (Ss)^ ensures that ()2.93p 
does not ask | sinc(jt| to be larger than 1 for any (S) that satisfies (|2.75p . so all the ellipses 
of (j2.77p are included. Every point on the boundary of the set of (S) described by (I2.75p 
is on the boundary of a positivity domain. This completes our proof that the compatibility 
domain and the intersection of the positivity domains both are the set of (S) described by 
the inequality ()2.75|) . 
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III. GENERAL FORMS 

Consider a quantum system described by x matrices. The NxN Hermitian matrices 
form a real linear space of A^^ dimensions with inner product 

N 

{A, B) = Tt[A^B]= Y: Al^B,,. (3.1) 

j,k=i 

Taking A^^ linearly independent Hermitian matrices that include the unit matrix 1, orthogo- 
nalizing them with a Gram-Schmidt process using the inner product 1)3.111 . starting with the 
unit matrix, and multiplying by positive numbers for normalization, yields A^^ Hermitian 
matrices F^q for /x = 0, 1, . . . A^^ — 1 such that Fqq is 1 and 

Tr[F,,oi^.o] = iVV- (3-2) 

Every NxN matrix is a linear combination of the matrices F^q. 
A state of this quantum system is described by a density matrix 

A^ 



P = 4 I 1 + E /^^^o I . (3.3) 



Equations ()3.2p imply that 
for /i = 1,2,. . . A^2 _ gQ 



(F^o) = Tr[F^op] = (3.4) 



P = ^ (^1 + E {Fao)Fao j ■ (3.5) 

Knowing p is equivalent to knowing the A^^ — 1 mean values (-F)io) for = 1,2,... A^^ — 1. 
The state is described either by the density matrix or by these mean values. We can see 
how the state changes in time by learning how these mean values change in time. 

Suppose this first system is entangled with and interacting with a second system described 
by M X M matrices. Let Fq^ for /i = 0, 1, . . . — 1 be Hermitian M x M matrices such 
that Fqq is 1 and 

Tr[Fo^Fo,] = M5^,. (3.6) 

The combined system is described by NM x NM matrices. Every NM x NM matrix is a 
linear combination of the matrices F^o ® Fo,y which are Hermitian and linearly independent. 
We use notation that identifies F^o with F^o ® 1 and Fq^, with 1 (g) Fqi^ and let 

F^u = F^o ® Fo,. (3.7) 
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For these NM x NM matrices 

1i[F^uFa,p]=NM 5^^5,p. (3.8) 

In the Heisenberg picture, evolution produced by a Hamiltonian H for the combined 
system changes each matrix F^^, to a matrix 

e'^*F^.e~*^* = E E (3-9) 

a=0 /3=0 



with real t^u;ap- Since 



Tr 



the t^y^ap form an orthogonal matrix, so t'^p.^y is t/^^ja/? and 

Af2-1 A/2-1 

e-^^*F«^e^^* = E E (3-11) 

Since Fqq is 1, 

Forming an orthogonal matrix is not the only property the t^y^ap need to have. They must 
also yield 

e'^*F^,e-*^* = e'^*F^oe-'^*e'^*Fo.e-'^* (3.13) 

and the same with t changed to —t. 

The mean values {F^q) for yU = 1, 2, . . . A^^ — 1 that describe the state of the first system 
at time zero are changed to the mean values 

(F^o)' = (e^^*F^oe-*^*) = + E ^MO;«o(i^ao) (3.14) 
that describe the state of the first system at time t, with 

= E E ^^^0■,»^B{Fa/3)■ (3.15) 
a=0 13=1 

Mean values (Fao) that describe the state of the first system are in equation ()3.14j) but not 
in ()3.15j] . We consider the ci^, as well as the tfj^o-ao to be parameters that describe the effect 
on the first system of the dynamics of the combined system that drives the evolution of the 
first system, not part of the description of the initial state of the first system. 



Tr[F^,F„^], (3.10) 
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The density matrix p of equation ()3.5|) that describes the state of the first system at time 
zero is changed to the density matrix 

P' = ^ (^1 + E^i^MO^o j (3.16) 
that describes the state at time t. Equations (|3.14|) imply it is 

P' = 1 + E d.F^o + E (^-o) E U,aoF,o . (3.17) 

\ ^l=l a=l fj,=l J 

Equation (|3.17|) for p' can be obtained another way. In the Schrodinger picture the density 
matrix 

n = ^ 1 + E iFc.o)F.o + E E {Fap)Fc.p (3.18) 

^^^^^ y a=l a=0 /3=1 / 

that represents the state of the combined system at time zero is changed at time t to 



e-^'^^ne^^* = ^ 1^1 + E (^-o) E Uc.oF,o 

+ E (^"o) E E 

a=l fi=0 u=l 

N^-lM'^-l Af2-1 

+ E E (^°/3) E 
a=0 13=1 fJ.=l 

+ E E E E (3.19) 

0=0 13=1 fi=0 u=l J 

according to equations ()3.11|) . Taking the partial trace of this over the states of the second 
system eliminates the F^y for v not zero and gives equation (j3.17p for the density matrix of 
the first system at time t with equations (j3.15|) for the d^. Since this involves working with 
the larger system longer, it does not appear to be the easier way to actually do a calculation. 

The map from density matrices ()3.5|) at time zero to density matrices ()3.17|) at time t 
holds for all the varying mean values (-Fao) that are compatible with fixed mean values {Fap) 
in the in describing a possible initial state for the combined system. We will refer to them 
as compatible (-Fqo)- Almost all initial states of the combined system allow the compatible 
{Fao) to vary as A^^ — 1 independent variables. We will consider only those initial states. 

The map of density matrices extends to a linear map of all x matrices to N x N 
matrices defined by 

1' = 1 + E ^f^F^o , F^Q = E itJi.O;aoF^Q. (3.20) 
11=1 ^=1 
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It takes the density matrix ()3.5p to the density matrix ()3.17p for each of the varying com- 
patible (Fao). It takes every Hermitian matrix to a Hermitian matrix. 

The latter property alone is the foundation for basic forms of the map. This statement 
is independent of our other considerations. 

Lemma: If a linear map Q ^ Q' of N x N matrices to N x N matrices maps every 
Hermitian matrix to a Hermitian matrix, then in the description of the map by 

N 

Q'rs — X! ^rj;skQjk (3-21) 
j,k=l 

the iV^ X A^^ matrix B is uniquely determined by the map and is Hermitian, 

Kj;sk = I3sk;rj, (3.22) 

and there are N x N matrices C{n) for n = 1, . . . A^^ such that 

Q' = Y,Cin)QC{ny - E C{n)QC{ny (3.23) 

n=l n=p+l 

for all Q, and 

Ti[C{myC{n)] = (3.24) 

for m ^ n, for m,n = 1, . . . N"^. 

Proof: Let Ejk be the N x N matrices defined by 

[Ejk]lm = ^Ij^mk- (3.25) 

Clearly Ej,^ = Ekj. If the map takes every Hermitian matrix to a Hermitian matrix, then 
{R.e[Ejk]y and (Im[£'jfc])' are Hermitian and 

{{E^kYV = {{MEjk])' + ^{MEJk]yV = iMEjk])' - ^{ME,k]y = {E]J. (3.26) 
Equations (HT^ and (IT^ give 

[-Sj^Jr-s — ^ ^ Brl-^sm^lj^mk — Brj.sk (3.27) 
l,m 

which shows that the map determines a unique B, and with 

(E;,)t = ^Ejj = E',^ (3.28) 
implies that B*j.j.f. = Brk;sj which is the same as equation ()3.22|) . 
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Since B is Hermitian, it has a spectral decomposition 



Af2 

B =Y1 \n\n){n\ 



(3.29) 



n=l 



where the \n) are orthonormal eigenvectors of B and the A„ are eigenvalues. The A„ are real, 
but they are not necessarily all different, non-zero, or non- negative. We label them so that 



An > for n = 1, . . .p 



A„<0 for n = p+l,...N^. (3.30) 



Then 



Af2 



Brj;sk = J2y\'^n\{rj\n){n\sk)^J\Xn\- J2 y\^n\{rj\n){n\sk)^J\Xn\. 

n=p+l 



n=l 



Let 



C{n)rj = ^/\K]{rj\n). 



Then equation ()3.2H) is 



Ar2 



n=l jk n=p+l jk 

SO the map is described by equation ()3.23|) . and 



(3.31) 



(3.32) 



(3.33) 



Tr[C'(m)tC(n)] =^C'(m):^.C'(n),, = 5:.\/|A^|(m|rj)(rj|n)v'|An| = \Xn\Hn) (3.34) 

rj rj 

which is zero for m 7^ in accord with equation ()3.24p . 
This completes the proof of the Lemma. 

The maps we are considering, those described by equations ()3.2()j) . have the additional 
property that 

Tr g' = Tr g (3.35) 

for every Q. This implies that 

Ar2 



P 

1 

n=l 



J2C{nyC{n)- C{n)^C{n) = l 

n=p+l 



(3.36) 



Tr g = Tr g' = Tr 



because 

YCinyC{n)- Y C(n)tC(n)|g (3.37) 

\,n=l n=p+l 

implies that in the linear space of NxN matrices with the inner product defined by the trace 
as in ()3.1|) . the difference between the two sides of equation ()3.36|) has zero inner product 
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with every matrix Q and therefore must be zero. From equations ()3.25p and ()3.27|) we see 
also that the trace-preserving property described by ()3.35|) imphes that 

^ Brr,rk = Tr = Tr [E.^] = 5,^- (3.38) 

r 

Conversely, either equation ()3.36|) or ()3.38|1 imphes that Tr Q' equals Tr Q for every matrix 
Q. From equation ()3.38|1 we see in particular that Tr B is A^. 

IV. DISCUSSION 

In the light of understanding gained here, it is easy to see the errors in arguments that 
a map describing evolution of an open quantum system has to he completely positive. One 
argument uses the fact that a map for a system A is completely positive if and only if it 
is the contraction to A of unitary evolution of a larger system AB in which A is combined 
with another system B and the density matrix for the initial state of AB is a product of 
density matrices for A and B. That is clearly not necessary. 

Another argument uses the fact that a map for a system A is completely positive if and 
only if the product of that map with the identity map for another system C yields a map 
for the combined system AC that takes every positive matrix for AC to a positive matrix. 
The argument says this is the way to satisfy the physically reasonable requirement that the 
description of the evolution of A must allow A to be accompanied by another system C that 
could be entangled with A but does not respond to the dynamics that drives the evolution 
of A. If the map for A is a contraction to A of either unitary evolution or a completely 
positive map for a larger system AB in which A is combined with another system 5, then 
the evolution of B is generally not described by the identity map, so C is not B. The 
accompanying system C must be a third system. The physically reasonable requirement 
can be satisfied very simply for the kind of maps we have considered. If the map for AB is 
completely positive, its product with the identity map for C yields a map for the combined 
system ABC that takes every positive matrix for ABC to a positive matrix. 

Mathematically, a map of states for a subsystem A can be constructed from (1) a map 
that takes density matrices for A to density matrices for the entire system AB at time zero, 
followed by (2) unitary Hamiltonian evolution from time zero to time t for AB., and finally 
(3) the trace over the states of B that yields the density matrix for A at time t. The broad 
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class of maps obtained this way is known to include maps that are not completely positive 
and in fact maps that do not take every positive matrix to a positive matrix. That all 
depends on the first step, the map that assigns density matrices pab to density matrices pA 
at time zero. Pechukas has shown that if A is a qubit, the only linear assignment of 
density matrices pab that applies to all density matrices pA, and gives back unchanged pA 
in the trace over B at time zero, is 

PA^ Pa® Pb (4.1) 

with Pb fixed. We prove this for any quantum system in an Appendix. Pechukas concludes 
that in general, when product assignments ()4.1|) do not apply, maps have to act on limited 
domains. This does not depend on the unitary evolution of AB from time zero to time 
t. When a product assignment ()4.1|1 is the first step, the map made in three steps is 
completely positive; if a map made this way is not completely positive, its domain must 
be limited. There has been debate whether any except the completely positive maps can 
describe physical evolution 15, IgI ]. 

Which do describe physical evolution? What is needed for one of these maps to describe 
evolution of states of A caused by dynamics of ABl If the map is meant to apply to a 
set of pA that all evolve in time as a result of the same cause, the pab assigned to these 
Pa should not differ in ways that would change the cause of evolution of the pA- If they 
did, we would say that different pA are being handled differently and that their evolution 



should be described by different maps. Pechkas jlj] considers the case where A and B are 
qubits and a product pab is assigned, as in (jUTj), to each of four selected pa, with a different 
Pb for each of the four pA- This yields a map that takes every mixture of the four pA to a 
density matrix. Pechukas observes that the large set of maps obtained this way must include 
many that are not completely positive and many that take density matrices outside the set 
of mixtures to matrices that are not positive. However, the pab assigned to each different 
mixture generally gives a different density matrix for B in the trace over the states of A. 
Each different state of A is coupled with a different state of B. Does this mean it is handled 
differently? If a map is meant to describe evolution that has a definite physical cause, does 
Pechukas have a single map that acts on a set of states; or a set of maps, each acting on a 
single state? 

In the compatibility domain that we describe, the evolution of all the states is clearly the 
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result of the same cause. It can be described by a single map that has physical meaning. 
Working with mean values helps make this clear. We do not need a complete description of 
the state of AB at time zero. It does not need to stand alone, independent of the unitary 
evolution, and accommodate any unitary evolution. The compatibility domain depends 
on the unitary evolution. In our example, the compatibility domain depends on the mean 
values that are the parameters Oi and 02. That these mean values are the relevant parameters 
depends on our choice of Hamiltonian. The compatibility domain is unlimited when oi and 
02 are zero. Then the map is completely positive, but that does not require an initial state 
described by a density matrix that is a product. 



APPENDIX: GENERALIZATION OF PECHUKAS' RESULT 

Theorem. If a linear map applies to all density matrices pA for a subsystem A and assigns 
each Pa a density matrix pab{pa) for the combined system AB so that 

Trij[pAB(pA)] = Pa, (A.l) 

then, for every p^, 

Pab{pa) = Pa® Pb (A. 2) 

with Pb a density matrix for the subsystem B that is the same for all pA- 
Proof. The first step, which Pechukas jl^ did, is to show that every pure-state density 
matrix pA is assigned a product density matrix, as in ()A.2|) . with ps possibly different for 
different pA- For completeness we include a slightly different presentation of this step. If 
Pa represents a pure state, there is an orthonormal basis of state vectors \ipj) for A, with 
j = 1, 2, . . ., such that pA is \ipi){'ipi\. We combine these with orthonormal state vectors 
for B to make an orthonormal basis of product vectors \ipj(t)k) for AB. Since Pab{pa) is 
positive, each {4'j(pk\PAB{PA)\'^j(pk) is non-negative and, from ()A.1|) . if j is not 1, 

(^j0fc|PAB(pA)|V^j0fc) < {lpj\TTB[pAB{pA)]\'ipj) = {ipjllpl) = 0. (A.3) 

Since Pab{pa) is positive, it is the square of a Hermitian operator. Thus we see that 
PAB{PA)\'^j<Pk) is zero if j is not 1 and 

{tjjj(j)r\pAB{pA)\^k(f)s) = Sji5kl{tpl(l)r\pAB{pA)\tpl(j)s)- (A.4) 
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Let 

PsiPA) = T^rA[pABiPA)]- (A. 5) 

Then 

Pb{pa) = (^iIpab(pa)I^i)- (A.6) 

and 

Pab{pa) = ®Pb{pa)- (A.7) 

That completes the first step of the proof. 



The second step, which completes the proof of the t 



leorem, is to show that is the 



same for all pure-state density matrices pa- Pechukas jl4| did this for the case where A is a 
qubit. We show that the proof an be easily extended to any quantum system. Suppose \ipi) 
and \ip2) are orthonormal state vectors for A. Let 

{ip^) = cosa It/'i) + sinae*^|'?/'2) 

\iPg) = sina IV"!) - cosae*^|?/'2). (A. 8) 

Then l^ps) and \ip4) are orthogonal, lip^) and l^pe) are orthogonal, and |('?/'i|'?/'3)p, |('?/'i|'?/'4)p, 
liHi^-s)?, \{Hi^5)?, KHA)?, mm?, and \{tP,\A)? are all 1/2. The length 
of each vector {ipk) is 1, so \'ipk){4'k\ is a pure-state density matrix for A. The map assigns it 
a product density matrix 

PABi\i^k){i^k\) = \^k){^k\ ® pBik) (A.9) 

as in ()A.7|1 with pB^k) short notation for PB{\4'k){4'k\)- 
Since the map is linear, it follows from 

^miH + \i^2){H) = li\i^3){H + |^4)(^4|) (A.io) 

that 

Pi?(l) + \^2){H Pb(2)) = limiH Pb(3) + \^,){H Pb(4)). (A.ll) 

Taking partial mean values {ipi\ ■ ■ ■ I'ipi), (V^2| • • • \i'2), (V'sl • • • iV^s) of this last equation ()A.11|) 
yields three equations that imply Pb(1), P_b(2), p_b(3) and Pb(4) all are the same. Doing 
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everything starting from ()A.10|) again with 1,2,3,4 changed to 3,4,5,6 shows that Pb(3), 
Pb(4), Pb(5) and Pb(6) all are the same. Any state vector for A is in a subspace spanned 
by \ipi) and a vector |?/'2) orthogonal to so \ip5) with fixed \ipi) and varying a, [3 and 
can represent any pure state for A. If p^ represents a pure state, Pb{pa) is the same 
as p_b(1), so ()A.2|1 holds, with the same pe, for all pure states of A and, therefore, for all 
mixtures as well. This completes the proof of the theorem. 
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